hep-ph/Ol 10331 



TTPOl-23 
NIKHEF 01-015 
October 2001 



Next-to-Next-to-Leading Order QCD Corrections 

^ ; to the Photon's Parton Structure 

o ' 
o 

(N 
> 

o 

■ S. Moch", J. A.M. Vermaseren^ and A. Vogt^ 

(N 

. "'Institut fiir Theoretische Teilchenphysik 

Universitdt Karlsruhe, D-76128 Karlsruhe, Germany 

g ; ^NIKHEF Theory Group 

Kruislaan 409, 1098 S J Amsterdam, The Netherlands 



ix. 



X 



Abstract 

The next-to-next-to-leading order (NNLO) corrections in massless perturbative QCD are derived 
for the parton distributions of the photon and the deep-inelastic structure functions and F2 ■ 
We present the full photonic coefficient functions at order aas and calculate the first six even-integer 
moments of the corresponding 0{aa^) photon-quark and photon-gluon splitting functions together 
with the moments of the aa^ coefficient functions which enter only beyond NNLO. These results 
are employed to construct parametrizations of the splitting functions which prove to be sufficiently 
accurate at least for momentum fractions x ^ 0.05. We also present explicit expressions for the 
transformation from the MS to the DIS^ factorization scheme and write down the solution of the 
evolution equations. The numerical impact of the NNLO corrections is discussed in both schemes. 



1 Introduction 



The hadronic structure of the photon, in particular the deep- inelastic structure function F2 , has 
attracted interest since the early days of perturbative QCD. Indeed, the leading-order (LO) correc- 
tions to the 'pointlike' parton-model result were derived twenty-five years ago in ref. Q, and the 
next-to-leading order (NLO) contributions followed a few years later Q. These results, obtained 
in the framework of the operator product expansion (OPE) Q, were recast in the language of 
evolution equations for the photon's quark and gluon momentum distributions in refs. and Q. 
An error in the NLO photon-gluon anomalous dimension was corrected ten years ago ||7|, ^. Unlike 
the case of lepton-nucleon deep-inelastic scattering (DIS) [9-15], the next-to-next-to-leading order 
(NNLO) QCD corrections have not been addressed for the photon structure up to now. 

So far the measurements of F2 have been performed using the process e~^e~ — >e"'"e~ + hadrons 
at electron-positron colliders; see refs. for recent overviews. While data from LEP have greatly 
improved the situation, an accuracy comparable to that achieved in lepton-hadron DIS can only be 
envisaged if 67 collisions will be realized via laser back-scattering |17] of one of the electron beams 



of a future linear collider [^, |19|]. Important information on the photon structure can also be 
expected from photoproduction of jets at HERA, which has been treated at NLO thus far [20 - 23] . 
The extension to NNLO is under way, e.g., the two-loop matrix element required also for hadronic 



collisions have been derived in refs. [24 1 using the pioneering results [^5[ for the scalar double box 



diagrams; the corresponding results with one external photon will be available soon [p 



In this article we present, with one qualification concerning the splitting functions, the NNLO 
corrections for electron-photon DIS and the evolution of the parton densities of the photon in 
massless perturbative QCD. In Sect. 2 we extend the OPE analysis of the photon structure [§, ^ to 
the required accuracy: The partonic forward amplitudes for the scattering of a virtual photon (and 
of a fictitious scalar directly coupled only to gluons) off a real photon are expressed in terms of the 
anomalous dimensions and coefficient functions up to second order in the strong coupling constant. 
Our calculation of these amplitudes for the lowest six even-integer values of the Mellin variable, 
N = 2, . . . , 12, then facilitates the extraction of the corresponding anomalous dimensions (up 
to NNLO) and coefficient functions (up to the next-to-next-to-next-to-leading order, N'^LO). The 
results for these quantities in the MS scheme are presented in Sect. 3 in numerical form, together 



with a brief discussion of the actual computation which closely followed the lines of refs. [11, 12 1. 
The analytic expressions for these results can be found in Appendix A. 

In Sect. 4 we switch to the parton language and specify the dependence of the photon-parton 
splitting functions and the photonic coefficient functions on the renormalization and factorization 
scales. After recalling the general factorization-scheme transformation, we then derive the NNLO 
corrections in the DIS^ scheme and discuss the 'physical' kernel for the non-singlet evolution 
of the structure functions at large values of the Bjorken variable x. The NNLO solution of the 
evolution equations is also given in this section. Explicit x-space expressions for the photonic 
coefficient functions and the photon-parton splitting functions up to NNLO are presented in Sect. 5. 
For the splitting functions we have to rely on our finite- results of Sect. 3, thus we can only 
provide approximations analogous to those derived in refs. 0, 15[ for the three-loop QCD 



splitting functions. For the NNLO coefficient functions and the corresponding transformation 
of the splitting functions to the DIS^ scheme we present, besides the exact results (deferred to 
Appendix B for the latter quantities), also compact approximate expressions. In Sect. 6 we finally 
illustrate the numerical effect of the NNLO corrections on the evolution kernels and on the solution 
of the evolution equations. Our conclusions are presented in Sect. 7. 
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2 Moments: formalism and method 



The subject of our calculation is inclusive hadron production in unpolarized electromagnetic (e.m.) 
deep-inelastic electron-photon scattering, 

e(k)+j{p) ^ e{k') + X , (2.1) 

where 'X' stands for all hadronic states allowed by quantum number conservation. The hadronic 
part of the corresponding amplitude is given by the (spin-averaged) tensor 



{P,q) = ^Jd^ze'^-'{T,p\J,{z)MO)\T,p) 



1 ^0. , 1 



with 



, .2x , , 

d,,u = -9tiu -Pfj,Pu-^ - {PfiQu +PuqtM) ■ (2.3) 

Here \T,p) denotes the physical photon state (including a non-perturbative hadronic component) 
with momentum p, and represents the e.m. quark current, q = k — k' is the momentum 
transferred by the electron, = —q^, and x = Q'^ /{2p ■ q) is the Bjorken variable (0 < x < 1). 
The longitudinal structure function is related to the structure function Fi by Fl = F2 — 2xFi. 



The optical theorem relates the tensor W^^, in Eq. (2^) to the forward amplitude T^,y for the 
scattering of a virtual photon off a real photon, 

VJ{p,q) = iJd'ze''^^r,p\T{Mz)MO))\T^,P) ■ (2-4) 

This quantity represents a convenient starting point for practical calculations, due to the presence 
of the time-ordered product of currents to which standard perturbation theory applies. In fact, 
the operator product expansion for this product and the subsequent application of a dispersion 
relation to Eq. (^^) very closely follow the procedure for standard lepton-hadron deep-inelastic 



scattering discussed, for example, in refs. |Tn, O, 27] to which we refer the reader for details. In the 



leading-twist sector addressed in this article, the only, but crucial difference to the lepton-hadron 
case is the presence of the spin-A^ twist-2 photon operators 0, ^ 

q{Mi,--.Mjv} ^ pH^^l £)l^l2 . . . £)f^N-l pt^N^^ (2.5) 

and their coefficient functions C/^ in addition to the usual quark (flavour non-singlet and singlet) 
and gluon operators, Onsi Oq and Og, and their respective coefficient functions. D^' in Eq. ( p.5| ) 
denotes the covariant derivative, and F^^ represents the e.m. field strength tensor. The spin- 
averaged matrix elements of these (renormalized) operators are given by 

(r,pI0i^--'^-^ lr,p) = p^^^K^p/'n} A''^J^?) , a = ns, q, g, 7 , (2.6) 



where /i stands for the renormalization scale. It is understood in Eqs. ( p.5[ ) and ( |2.6| ) that the 
symmetric and traceless part is taken with respect to the indices in curved brackets. 
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Following the procedure for lepton-hadron DIS |10, 11, 27], the even- integer Mellin-A^ moments 
of the structm'e functions and in Eq. (2.2) 



Ti{x) = -Fi{x) 



can then be expressed in terms of the parameters of the OPE, 



2,L 



(2.7) 



(21 



a=ns,q,g,7 

Here and throughout the whole article we use the notation 



4tt 



Or 



a 
4tt 



(2. 



for the strong and electromagnetic coupling constants. The present study addresses the higher- 
order QCD corrections to the photon structure functions at the leading order of QED, al^. 
Consequently the quantities entering the r.h.s. of Eq. ( ^ ) are only needed at their respective lowest 
e.m. orders, i.e., Ar,y and Cj^p (p = ns,q, g) at a^^, and ^r,p and Cj_^ at a^^. 

The operators Oa in Eq. ( |2.6| ) mix under renormalization. Expressing the renormalized opera- 
tors in terms of their bare counterparts, this mixing can be written as 



(2.10) 



Here and in the next two equations the summation convention is used, and the range of all indices 
is as specified in Eq. ( |2.6| ) above. The anomalous dimensions jap governing the scale dependence 
of the operators Oa, 

-^Oa = -JafsOp , (2.11) 



din fi^ 

are connected to the mixing matrix Zap in Eq. ( |2.10| ) by 

d 



z 

din /i^ 



^ {Z ^)a'l3 ■ 



(2.12) 



Keeping only those terms which are relevant for the structure functions i^g^ at order a],^, the 
matrices Z and 7 take the form 



Z 



^ Zns 








y 







Zqq 


Zqg 


Zoci 







Zgq 




Zg.y 




V 








1 





(2.13) 



and 



fin 



7, 



(0) 



with the perturbative expansions (where kj 

00 

7p(p') = E ^Sv') ' 









k 


\ 


7qq 


7qg 


kq 




7gq 


7gg 


kg 













/ 


0) 









(2.14) 



00 

1=0 



^em ki^ ^ 



(2.15) 
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In order to make practical use of Eq. ( 2.12 ) a regularization procedure and a renormalization 
scheme need to be selected. We choose dimensional regularization and the modified |^] minimal 
subtraction scheme, MS — the standard choice for modern multi-loop calculations in QCD. 
For this choice the running couplings in D = 4 — 2e dimensions evolve according to 



daj 
din ^2 



-e Oj + I3j{aj 



3 



em, s 



(2.16) 



where /?em and (3^ denote the usual four-dimensional beta functions of QED and QCD, respectively. 
/3em does actually not enter the present calculation; for [5^ we employ the standard notation 



/3s 



-/3o al - I3ial- /?2 a 



(2.17) 



with (3q = 11 — 2/3 N{, where A'^f stands for the number of effectively massless quark flavours. 

Inserting Eqs. ( |2.13| )-( ^.17| ) into Eq. ( |2.12D and solving for the photon-quark and photon-gluon 
renormalization factors Z^s^, Z^^ and Z„^ we obtain, up to the desired order a, 



a2 

""erne's ' 



1 



1 



1a:(o)^(o) 

n ^ns /ns 



+ < 



1 /it.(o) f^w^i^ - ^R^m^i'^) 

g3 I g ns \^ /ns j gMO^'ns ms 



g2 I 3 'i'^ ''i'^ g ns /ns gMU"'ns 



3e 



(2.18) 



-1 7 



■ fef + as 



+ < 



2 q 'qq 

/qq 



+ — fcW 
+ 2e^^ 



_L i 1.(0)^(0) (0) _ i«^i,(o) (0) 

^ g'^q 'qg 'gq g'^Oft'q 7qq 

+ g2 ^qq + gN ^qq + g'*^g ^qg g/^o^q 



3e ^ 



(2.19) 



^em ^g7 



+ Os 



+ < 



1 n 



A;(0)^(0) 

2 q 'gq 



+ 2e s 



.3 \gN ^gq V^qq + ^gg ) g^^ON ^gq 
+ g2 ^gq + gN ^gq + g'*^g ^gg g/'oKg 



+ -4^) 

3e s 



(2.20) 



The photon-parton anomalous dimensions /cnsj ^q and /cg can thus be read off order- by-order from 
the terms of the corresponding renormalization factors, while the higher poles in 1/e in Eqs. 
( 2.18 )-( 2.20 ) can serve as checks for the calculation. The coefficient functions in Eq. (p.^), on the 
other hand, have an expansion in positive powers of e, viz 



Ci,p 



*i7 



k2 (1 - <5, 

00 



00 

pgv + E "-s ("^1! 



(0 
>p 



J-1 



,(0 

'i,7 



(0 



1,7 



(2.21) 



where i = 2, L and p = ns, q, g. Also here only the lowest-order terms in Oem have been retained as 
discussed below Eq. (|2.8| ), thus, like the anomalous dimensions 7p(p') in Eq. ( 2.15| ), the coefficient 
functions Q^p in Eq. ( p.21| ) are just the standard QCD quantities entering lepton-hadron DIS. 
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Due to the partly non-perturbative physical photon state \T,p), Eqs. (0) and (^) are not 
accessible to a perturbative computation. However, as the OPE represents an operator relation, 
the anomalous dimensions ( p.l5 ) and the coefficient functions (2.21) do not depend on this state. 
Hence, again closely following the procedure []lO| , |ll], 12| for lepton-nucleon DIS, the calculation can 
be performed using a partonic photon state \ j,p)- Instead of Eq. (|2.4|) we thus consider 



(2.22) 



VJ{p,q) = z d^ze^'^^^,p\T{J^iz)MO))\^,p) 



At leading- twist accuracy the decomposition of T^J into T2^j and Tl^^ analogous to Eq. (2^) is 
provided by 



/ 3 - 2e q 



p^p'' + 



2e {p^ 



2e 



(2.23) 



The A^**^ moments are obtained from Eqs. ( 2.23 ) by applying the projection operator Vn |31] 



^ i,7 



, flg , Clem ) £ ) 



(2.24) 



qlt^i . . . gMJv} 
iV! 



^p^^^ . . .dp' 



Ti^^{x,Q , Cls, Ctemj c) 



p=0 



where g'f'^^ • • • q^^'^^ is the harmonic, i.e., the symmetric and traceless part of the tensor 

This operator does not act on the coefficient functions Ci^a and the renormalization constants 



Zap in Eq. (2.10), which are functions only of A^, as, a and e. It does act, however, on the bare 
matrix elements ^ (defined analogously to Eq. ( p.6D ) and eliminates all diagrams containing 
loops, as the nullification of p transform these diagrams to massless tadpole diagrams which are 
put to zero in dimensional regularization. This removes the operator matrix elements A^^-, P = ns, 
q, g, which only start at the one- loop level. Hence only the matrix elements A^^^'^'^'^ of the photon 
operators (|2.5|) remain. These matrix elements are given by 



4Y'-^^(e) = (1 - e) • const^ , 

where the factor (1 — e) arises from the number of photon polarizations va. D = A 
We thus arrive at 



(2.25) 
2e dimensions. 



rpN ( Q 



E 



a=ns,q,g,7 



a7 



(2.26) 



with i = 2, L. This relation, after expanding in powers of aem, and e, provides a system of 
coupled equations which can be solved for the photon-parton anomalous dimensions and the photon 
coefficient functions. In particular, by computing T^t^ to the order Oem > we can derive the desired 



coefficients kns , kq , and c-^ in Eqs. ( 2.15| ) and ( 2.21 ). The quantities kg , on the other 
hand, cannot be determined in this manner since the gluonic coefficient function in Eq. ( 2.21| ), 
unlike its quark counterparts, starts at order Os only. This problem is overcome by considering, in 
addition to Eq. ( 2.22| ), another unphysical Green function T'^''^ where the virtual-photon probe is 
replaced by an external scalar field cp coupling directly only to gluons, see below. 
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The expansion of Eq. ( ^.26 ) to order can be written as 



~ ^ CLcmO-s^e^ ) ' (2.27) 

The factor Se = exp [e{ln(47r — 7^}], where 7^^ denotes the Euler-Mascheroni constant, is an artefact 
of dimensional regularization kept out of the coefficient functions and anomalous dimensions in the 
MS scheme. The expansion coefficients Tf}^ can be decomposed into flavour non-singlet (ns) and 
singlet (s) pieces, 

= S^tI^"^ + ^) 5,tI^^ (2.28) 

which collect the contributions proportional to the respective combinations of quark charges, 

Sns = 3N, ((e^) - (e^)^) , 6^ ^ 6^ = SN^ (e^) ^ 3Efii e|, . (2.29) 

The pure-singlet (ps) contribution defined in the second line of Eq. ( |2.28| ) starts only at order a^. 
Consequently the anomalous dimensions k^} and k^^ in Eq. ( |2.15|) are identical for / < 1 except 
for their obvious charge factors. The same applies, for / < 2, to the non-singlet and singlet photon 
coefficient functions in 

^ = ^T + ie")^ ■ (2.30) 
The corresponding decomposition for the hadronic quantities reads 

(/) = (0 + (/) ^(0 = ^(0 + g(0 ,2 31) 

/qq /ns ^ /ps i H,q S,ns ^ «,ps ' V^-*-"^/ 

where the pure-singlet contributions are non-vanishing for / > 1 (2) in the first (second) relation. 



The explicit expressions for the first two singlet expansion coefficients in Eq. ( |2.28| ) in terms of 
the anomalous dimensions and coefficient functions are given by 

6.T^^'' = ^<)+c«'^ (2.32) 



and 



The corresponding non-singlet relations are obvious as discussed below Eq. ( p. 29 ), hence they are 
not written out for brevity. The contributions at the order read 

. y(2),ns _ 1 fl^(O) (0) (1)1 , 1 k(0) (2) , . iyt(0)^(0)a(l) 

^ '^L,7 ^ '^ns "L,ns ^ ^ "'^ ^'"^ 2 '^'^ -^'^ ' 
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r rp{2),ns 



1 /li.(o) - -BnUo) (0)' 

1 / 1^(0)^(0) Ji) + 1a;(0)^{i) + iA;(i)'YW - i/3nA:(i)l 
g2 I 2 "'^ 'qq 2,q 3 6 g^o'^'ns r 

^ ^ii.(2) , i.{o)„(2) , i> (0) (0) (1) 



, „(3),ns . (0) (2) k(0) (0). (1) 

''2,7 "-ns "'2,ns 2 '^'^ 2,q 2 "'^ '^'1^ 2,q 



(2.34) 



and 



r 7^(2),ps 



r T^(2),ps 



g2 1 2 -^'gj e I "1 -^'P^ 2 ^ -^'S 2 "1 ^^'^ 



+ c 



q i,ps 



L,g "T 2 ""q ')'gq "L.g ' 

1 /ii.(0)^(0)^(0)\ , 1 /ii.(0) (0) (1) , k(i)^(o) , i 1.(0) (1)1 
g3 1 g'^q 'qg 'gq J g2 1 2 6^ 3 J 



(2.35) 



Finally we derive the corresponding expressions for the unphysical flavour-singlet Green function 
T'^^{p, q) mentioned below Eq. ( 2.26 ). After application of the projection operator Vn in Eq- ( 2.24| ) 
the moments of T'l''^ can be written as 



-^^7 1 ^2 '«s,ae 



Q 



5^ ^^,a i 2") '^S) Q^eim £ j ■^(^ ('f^S) Oem) A 
=q,g,7 / ^ 



A'', tree 
77 



(2.36) 



The coefficient functions C^^a have an expansion analogous to Eq. (|]2l|), but with C^^g = 0(1) and 
C(f),q = 0{as). It is understood in Eq. ( p. 36 ) that the external gluon operator Gi^i^G'^'^ employed 
for the coupling to the scalar field (p has been renormalized according to [32| 



flU 



Zg2 (G^'^G^,)^^'^^ + 



1 - /?(Qs)/(eas 



(2.37) 



where the dots indicate mixing with unphysical operators which give vanishing contributions under 
the conditions of the present calculation. Expanding Eq. ( p.36| ) analogously to Eq. ( p. 27 ), the 
coefficients T^/l,, I = 1,2 read 



^(1) 



1 n 



(2.38) 



and 



.(2) 
<A.7 



e3 16 

1 n 



1 n 

e 



N 7gq l^Tqq + 7gg J gPoKq 7gq 

A;(o)^(o)cW + ifcW-yWcl^) + -k^o) (1) , 1^(1)^(0) , 1^(1)^(0) _ Ignk^^A 
N 'gq '-./..g ^ 2 'I 'qq 0,q ^ 3 q 'gq e q '^m 6 s j 

i.(2) , ^.(0) J2) , 1^.(1) Ji) , + iiu(o)^(o)„(i) 1,(0) (0) (1) 1 ^2 3q) 



In Eq. ( |2.39| ) we have not written out the term fixing the unphysical coefficient function 
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3 Moments: calculation and results 



The calculation of the moments ( |2.26| ) and ( 2.36| ) of the Green functions T^^ and T'^^ can be 



performed quite analogously to that of T'^s and T'*'^ in refs. [11, (where these quantities are 
denoted as T^"^^"^ and T^'t'^'t'')^ to which the reader is referred for a more detailed discussion. Indeed, 
the Feynman diagrams contributing to T'^'^ (^x'^"') derive from a subset of those for T'^s (T'^^): 
at three loops 117 (57) out of 366 (7162) diagrams contribute, respectively, using the counting of 



refs. 1 11, Despite this reduction we are not able to compute higher moments than obtained for 



the hadronic case in ref. [12|, as some of the most storage-consuming diagrams remain. 



The moments = 2, . . . , 8 have been calculated from scratch. The diagrams are generated 
using a special version of QGRAF |33|. The actual computation is done using optimized FORM 



1 34] programs using the MINCER package [^] for the scalar three- loop integrals. The moments 
iV = 2,4 of and iV = 2 of T'^'^ have been computed in an arbitrary covariant gauge, i.e., 
keeping the gauge parameter in the gluon propagator, i [—g^^ + (1 — O'i^'f / i'f + + ^ 

a free parameter. The explicit cancellation of the gauge dependence in the anomalous dimensions 
and coefficient functions provides an important check of the results. The direct calculation of the 
moments A^=10 and A^ = 12 is rather time-consuming and requires running FORM on a computer 



with 64-bit architecture instead on a standard PC |12]. Therefore we have for these moments (and 
as a further check also for the lower moments) made use of the diagram database of ref. |l^ by 
replacing the colour factors of T"^^ and T'^^ by those for our cases and then re- assembling the 
integrated results of all diagrams, thus sidestepping the involved parts of the computation. 

From the results of these calculations the 0{a^^a'^) photon-parton anomalous dimensions and 



photon coefficient functions at = 2, . . . , 12 can be obtained by means of Eqs. (2.27), (2.2. 



( 2.321) ^( 2.35 ) and ( |2.38| )-( 2.39| ). Here we present the results in numerical form; the full analytic 



expressions can be found in Appendix A. The non-singlet (p = ns) and singlet (p = s) photon-quark 
anomalous dimensions read 

(Sp Ocm)"^ fc^"^ = - 1.333333333 - as 7.308641975 

2 



+ at ( -3.368998628 /ps iVf - 86.97527479 + 1.470507545 N{ 

p "em j ^ 



{Sp flem)"^ C"'^ = - 0.7333333333 - Os 8.343259259 



+ ( -0.4223537997 /ps A^f - 102.8310333 + 1.477370645 A^f 



{Sp aem)~^ = - 0.5238095238 - as 8.525587590 



p 

,2 



+ < ( -0.1483655726 /ps A^f - 109.2777760 + 1.656530217 A^f 



{Spacmy^k^~^ = -0.4111111111 -as 8.343215919 



+ as^ ( -0.06961022336 /ps A^f - 111.1669036 + 1.695497529 A^f 

=10 _ 

p 



i^p fflcm)"^ k^~^° = - 0.3393939394 - as 8.045400445 



+ ( -0.03823222708 /ps A^f - 111.0346700 + 1.670607074 A^f ) 
(5pacm)~^ = - 0.2893772894 - as 7.720564967 (3.1) 

+ a^ ( -0.02324832852 /ps A^f - 109.9430144 + 1.619078961 A^f ) , 

where the factors 5ns and 6s have been defined in Eq. ( |2.29| ), and 
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The corresponding results for the photon-gluon anomalous dimensions are given by 
((5s Oem)"^ fe^"^ = as 1.975308642 + ( 31.41971923 + 5.157750343 Nf ) 
(4 aem)"^ =^ = as 0.8743703704 + ( 23.94271102 + 1.108863649 N{ ) 
((5s aem)~^ A;^=^ = as 0.4439549365 + ( 15.65166504 + 0.6959529917 iVf ) 
((5s aem)~^ k^^^ = as 0.2669916389 + ( 10.96606146 + 0.4981961912 iVf ) 
((5s ttem)"^ /cf = as 0.1780490673 + ( 8.160308004 + 0.3790604508 iVf ) (3.3) 
{5saem)~^k^='^^ = as 0.1271644566 + 0^(6.348294844 + 0.3002739601 ATf) . 
For the scale choice fi^ = the photon coefficient functions for F2 and Fl sX N = 2, . . . , 12 read 
('5paem)"^C2Pf=^ = -1 - as 5.54162471 

+ al ( -5.818386905 /ps iVf - 83.61538634 + 12.00984563 A^f ) 
(<5p aem)"^ 4;^^"^ = - 1.477777778 - as 12.20420558 

+ 4 ( -12.34431320 /ps iVf - 270.6890263 + 22.99272397 A^f ) 
{Sp aem)"^ 4;^^^ = - 1.410317460 - 16.88045086 

+ a^ ( -10.77343760 /ps iVf - 461.1411168 + 31.41856945 Nf ) 
((5p aem)"^ c^;^^^ = - 1.288174603 - as 19.73136088 

+ al ( -8.964615998 /ps Nf - 617.8190552 + 37.24672336 A^^f ) 

((5p aem)~^ cl'^^'^^ = - 1.172255892 - as 21.44335073 

+ al ( -7.480750665 /ps Nf - 741.1418081 + 41.20962391 A^^f ) 
{6pae^)-^cl'^=^^ = - 1.071686118 - as 22.46033065 (3.4) 

+ al ( -6.308511998 /ps Nf - 837.9184018 + 43.90927770 Nf ) 

and 

{Spaemy^cll^^^ = 1.333333333 - ag 12.54891027 

+ al ( 8.693227731 /ps A^f - 339.8638216 + 26.87433155 A^f ) 
(Spaeu,)'^ cll!^^^ = 0.5333333333 - as 6.551703704 

+ al ( -0.8773595064 /ps A^f - 208.9707023 + 16.63839609 A^f ) 
(Spaeui)'^ cl'H^^^ = 0.2857142857 - as 3.966771047 

+ al ( -1.452744800 /ps A^f - 138.7560249 + 10.96941052 A^f ) 
((5p aem)"^ cl]!^^^ = 0.1777777778 - as 2.670167875 

+ al ( -1.281927733 /ps Nf - 99.25459329 + 7.787570553 Nf ) 
((5p aem)"^ = 0.1212121212 - as 1.927770988 

+ al ( -1.056290000 /ps Nf - 74.93923167 + 5.838900780 Nf ) 
{Spa^^y^cl'f^^'^'^ = 0.08791208791 - as 1.461875421 (3.5) 

+ al ( -0.8667773331 /ps Nf - 58.86856783 + 4.557749952 A/^f ) . 
The additional terms for /x^ 7^ do not contain independent information, see the next section. 
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4 Parton distributions and evolution equations 



In this section we outhne the parton formulation of the photon structure (introduced in ref. ||5|) 
at the next-to-next-to-leading order (NNLO) of perturbative QCD. The number distributions 
of quarks and gluons in the fractional photon momentum x are denoted by qj{x,fi'j,fi'^) and 
g'^{x, /Xj, //^), where the subscript j indicates the quark flavour. |U/ represents the mass-factorization 
scale which, unlike in Sect. 2, is not generally identified with the coupling-constant renormalization 
scale denoted by fir from now on. The photon's quark and antiquark distributions are equal due 
to charge conjugation invariance. 

First we focus on the flavour-singlet distributions for which we use the notation 

) , s--E(^7 + ^7) = 2E^7' (4-1) 



where S"^ is the singlet quark density and A'f, as before, stands for the number of effectively massless 
quark flavours. As in some other equations below, the dependence on x and the scales has been 
suppressed in Eq. (|4.lD . At lowest order in the electromagnetic coupling Oem these distributions 
are subject to the evolution equations 

ASL = pf + P(^q'y , (4.2) 
a m fij 

where represents the Mellin convolution in the momentum variable, 

[a0b]{x) = f'^a{y)b(-) , (4.3) 
Jx y \yJ 

and the splitting-function matrices are given by 

p-f ^ ( 'Pn-y \ p ^ [ ^qq ^qg j _ (4 4) 

\ '^g7 / \ '^gq '^gg / 

The NNLO expansions of the photon-parton and parton-parton splitting functions read 

P^ (x, ttem, as, Luj = ttem P'^^ (x) + flemOs P^^^ (x) 

+ aemfls {P^^Hx) - 13oLrP^^\x)) +... 

p(x,as,L/j) = asP^^Hx) +al[p^'\x)-poLRP^^^{x)) (4.5) 
+ a! (P(2)(^) _ 2PoLrP^'Hx) - {i3iLr - P^L^P^^Hx)) +... . 



The notation for the coupling constants and /3o has been defined in Eqs. (|2.9| ) and ( p. 17 ) above, 
and we abbreviate the scale logarithms by 

Lm = , Lr = In^ . (4.6) 

The expansion coefficients in Eq. (4^) are related to the (even- integer A^) anomalous dimensions 
in Eq. {^J§ by 

fe?^ = - f^dxx^-'P^^ix) , 7^ = - [dxx^-'P^l^ix) . (4.7) 
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The x-space splitting functions, in turn, are uniquely fixed by the inverse Mellin transform if the 
anomalous dimensions are known for all even moments A^. 

The evolution equations for the non-singlet combinations of the photon's quark densities, viz 

= <ll - <l] ^ l<^/j<iVf (4.8) 

and linear combinations thereof, are obtained from Eqs. ( f4.2D and ( [4 .51) by replacing P and 
by their scalar non-singlet counterparts P^^ and P^y The hadronic quantities -P^'^ have been 
discussed, e.g., in refs. |13, 15]; the splitting functions P^^ are related to the non-singlet anomalous 
dimensions in Eq. ( |2.15| ) by 

3(e| - el)6-Mf = " [' dx x'^"' pg{x) . (4.9) 







Due to qj = there are no combinations evolving with the splitting functions Pj^ and P^s here, 
unlike for the nucleon structure. Hence the calculation of e.m. DIS in Sect. 2 and Sect. 3, if extended 
to all moments, is sufficient to fix all NNLO splitting functions relevant to the photon structure. 
Note also that differences of quark densities of equal charge, like d"' — s'^ , evolve like hadronic 
non-singlet distributions, i.e., without an inhomogeneous term in Eq. ( [4.2[) . 

The anomalous dimensions kp^^ , p = ns, q, g, for the next-to-leading order (NLO, I < 1) 
evolution have first been written down in ref. |^], their x-space counterparts Pp^\x) in ref. [^.Q 
Both results are correct except for the constant-TV (6(1 — x)) term erroneously included in kg'^^ 
{Pg^\x)), an error which has been corrected in refs. [^, The flem^s contributions Pp^^(x) to 
the photon-parton splitting functions are not completely fixed at present, as the moments > 14 
remain uncalculated. Following the lines of refs. [|^, |l^, [l^, we will employ our six moments ( p.l| ) 
and ( p^ ) to derive approximate results for Ppj (x) in Sect. 5 below. 

In terms of the parton distributions ( |4.1| ) and (O), the 0{al^) electromagnetic structure 
functions = ^ F-^ and ^2 = ^ F2 are given by 



«,7 



,7 r„2 „2> 



+ {e^) {ci^{x,acm,as,LM,LR) + d {x,as, Lm, Lr) ® q"'{fij, fif.) {x)j . (4.10) 
The first line in Eq. ( 4.10[) represents the non-singlet contribution involving the combination 

of quark densities. The x-space coefficient functions Ci^^{x) and Cj^p(x) are the inverse Mellin 
transforms of the e'' terms in Eq. ( |2.21| ) , and in the second line of Eq. ( [4.10| ) we have employed the 
matrix notation Cj = ( Cj^q, ) for the hadronic singlet coefficient functions. Up to order flem^s 
the expansion of the photonic coefficient functions takes the form 

^cTm Ci,'y{x, flcm, As, ^A/, Lr) = C^^^ (x, Lm) + AsCf,^^ (x, Lm) + O^C^f^ {x, Lm, Lr) 

E fcS;i(x) + E c^:f{x)L^\ . (4.12) 

1=1 \ m=l j 



^Oux notation for the anomalous dimensions fcp''^ differs from that for Kp^'" in ref. Q (where the analogue of 
Eq. (2.10) is written in terms of d/dlnfi instead of d/d\nfi^) by a factor -1/2. Our splitting functions Pp-.j{x) are 
larger by a factor 2'^^ than their counterparts fcp^(x) in refs. ^, 0, ^, where the expansion parameters are normalized 
as aj/2n, j = em, s, instead of Eq. (2.9). 
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The first-order contributions cj^^ liave been obtained in ref. Q by exploiting their simple relation 

to the NLO gluon coefficient functions c[^^ |2£]. Also the second-order coefficients c^^^ can readily 

be inferred from their gluonic counterparts c[^^ first calculated in ref. Q; explicit expressions will 
be presented in Sect. 5 below. The terms up to order Oemffls contribute to F{' and F.^ in the NNLO 
approximation. The fflem^s parts of Eq. (|4.12| ) only enter for dFi^2/dlnQ^ to this accuracy, hence 

(3) 

the scale-independent quantities c]^, i = 1,2, are not required. 

The coefficients cf'^\x) of the scale logarithms in the second line of Eq. ( |4.12 ) can be derived 
by identifying the results of the following two calculations of {d/d In Q'^Y Fi for / = 1, 2, 3 at 

= fjj = Q'^ : (a) with the scales identified in the beginning, using Eq. ( |2.16| ) for D = 4 and j = s 
together with Eq. ( [4. 21) , and (b) with the scales set equal only at the end, after differentiating the 
logarithms in Eq. ( ^.12| ). For the singlet contributions one finds 

,(3,1). ^ ^(0) ^ p(2) ^ ^(1) ^ ^ ^(2) ^ ^(0) _ ^(2),s 

cf^f'^ = ijcf ) ((P(0) - /?ol) P« + ® P(°)) + c« ® (P(0) - 2/3ol) Pf } 

cg3),. = icf)0P(o)^(p(o)-2/3ol)0Pf , (4.13) 

where cf^\x) = (C'^^^ (a;) , cf'gix)) = ((5(1 — x) , 0) is the parton model result for lepton-hadron DIS, 
and 1 represents the 2x2 unit matrix multiplied by (5(1 — x). The corresponding non-singlet results 
are obtained from Eqs. ( 4.13| ) by simply replacing all quantities on the right-hand-sides by their 
scalar non-singlet analogues. Finally the coefficients C^^ for fxj ^ \ir in Eq. ( ^.12| ) are obtained 
from these results by expanding as(M/) ™ terms of as{fi'^), 

eg (^. Lm, Lr) = eg (x, Lm, 0) - PqLr eg (x, Lm) • (4.14) 

The relations for the hadronic coefficient functions Ci{x) corresponding to Eqs. ( 4.12| )-( |4.14| ) are 
given in Eqs. (2.16)-(2.18) of ref. [|l4|. 

The decomposition ('factorization') of the r.h.s. of Eq. ( [4.10| ) into coefficient functions and 
parton densities is not unique beyond the leading order (LO), where the structure functions J-[^ , 
i = 1,2, are simply given by q^^'~' + (e^)S''''^^. Working, as always, at lowest order in Oem, 
the general factorization scheme transformation (corresponding to a finite renormalization of the 
operators in Eq. ( |2.6| ) ) of the singlet parton distributions ( [4.1D can be written as 

q"f = Z^ + Z0q^ . (4.15) 

Here Z'^ , like , is a two-component vector and a 2 x 2 matrix. It is sufficient to discuss the 
scheme transformations with all scales identified, hence we put /i^ = = in the following. By 
virtue of Eq. ( [4.2| ) , the evolution equations for the transformed distributions q^ then read 

dq^ /J.^ + Z^P^+f/Js^+Z^p) ^Z-^®(q^-Z^ 



dlnQ dtts \ da 

= P"' + P^q^ , (4.16) 
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where Ps = /3qcd has been defined in Eq. ( 2.17 ). The transformed coefficient functions read 



a 



1,7 



c, 



1,7 



Ci 



(4.17) 



The second term proportional to in Eq. (4.16) and the second relation in Eq. (4.17) are the 
standard expressions for scheme transformations of hadronic parton densities. 



Like their gluonic counterparts cfl, the photonic MS coefficient functions q-'i, are negative and 



(0 



■«,7 



singular for x ^ 1, the leading terms read 74/ln^'^^(l — x) with Ai > 0. Unlike the gluon case, 
however, the behaviour of Cj^-y enters the structure functions ( [4.1[1| ) directly, not via a regularizing 
convolution with a parton density. Hence, in the MS scheme, these singularities have to be com- 
pensated by the quark distributions which thus have to be rather different from the leading-order 
quantities at NLO and NNLO. Consequently the perturbative stability of the expansion in is 
more manifest in factorization schemes in which the photonic coefficient functions are removed 
as far as possible. The minimal modification of the MS scheme that removes (at fij = Q2) this 
function for the most important structure function F2 is the DIS^ scheme introduced in ref. Q 
(for a further discussion see also ref. Q]), where F2 is written as 



(4.18) 



with the standard MS hadronic coefficient functions C2,ns 3.nd C2. Eq. ( 4.1^ ) is obtained for the 
transformation 



^7,DIS^ 



''2,7 ~^ 







^DIS^ 



(4.19) 



Inserting Eq. ( 4.19| ) into Eq. ( |4.16| ) yields the NNLO singlet photon-parton splitting functions in 
the DIS-y scheme in terms of the MS splitting functions and coefficient functions]^ 



pDIS^ 

q.1 



,DIS^ 



Pr 



q7 



g7 



P. 



SI 



"s qq 



""S gq 



-2,7 
-2,7 



/p(l) 

\ qq 



.(l),s 



-2,7 



p(0) 

qq 



/?o 1) (S) (( 



,(2),s 



■2,7 



«'«4;i)} 

{pS S + 4? » (<=S' - <=S' » <=K) } ^ (4.20) 



The corresponding non-singlet relation is obtained from the first line of Eq. ( 4.20 ) by substituting 
the non-singlet counterparts of all quantities. Finally, due to Eq. ( [4.17 ), the photonic coefficient 
functions for J^^ in the DlS^y scheme coincide with those of —J'l in the MS scheme. 

The structure functions F,P'{x,Q'^) evolve approximately like non-singlet quantities at large x, 
since the photon's gluon density is much smaller than the quark distributions in this region and the 
pure-singlet contributions to both the splitting functions and the coefficient functions are strongly 
suppressed at large x (large A^), see, e.g., Eqs. ( [3.1D and ( ^^ ) for the photonic quantities. This non- 
singlet evolution can conveniently be recast in a form in which any dependence of the factorization 
scheme and the scale /x/ is explicitly eliminated, viz 



I,^{x,a,,Lii)+ /Ci(as,Li?)®^,l3(Q2) (x) . 



(4.21) 



The perturbative expansions of the 'physical' evolution kernels I^' and /Cj, (i = 1,2) for fif. = 



can be obtained by inserting Z'^ = d^-y and Z = Cj^ns into the non-singlet analogue of Eq. ( 4.161 ) 



^Note that a minus sign is missing on the r.h.s. of the second line of Eq. (3.1) in the journal version of ref. 
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yielding 



+ cgs ^ ^2} - ^ cQ- - A) cS^^ - {Pil^ - A) cS^,) cS J^} + . . . (4.22) 



and 



(4.23) 



The kernels and /Cj for ^ can be reconstructed from these relations by inserting the 
expansion of as{Q'^) as a power series in as(//^). Eq. ( [4.22| ) has been taken over from ref. the 
hadronic contribution (4.23) has been discussed up to order in ref. |39|. 

The LO splitting functions Pils7 and , and the combinations of splitting functions and 
coefficient functions enclosed in curved brackets in Eq. (4.22) and ( [4.23 ) represent factorization 
scheme invariants. Hence, together with the hadronic coefficient functions cj^g, the quantities 
-fns\ Cj-^'"^ and Pns-y have to be included for the NLO approximation of the photon's structure 
functions F^^^, i = 1,2. Correspondingly, the NNLO evolution additionally requires, besides q^g, 

the functions Pis\ cf'}^'^^ and Pils7- In general, the terms up to the orders a^^a^ {(^s^^) the 
first (second) part of Eq. ( ^ ) and a^^^a^"^ (a^ ) in the photonic (hadronic) coefficient functions, 
respectively, contribute to the N'^LO approximation together with the coefficients Pi<k in Eq. ( 2.17| ). 

We now turn to the NNLO solution of the evolution equation ( [4.21) . Here we confine ourselves 
to the (sufficiently general) case of a constant ratio expL/j = fi'j/fi'^. In this case the '[/-matrix' 
technique, developed for the NLO hadronic evolution in ref. and generalized to higher orders 



in refs. |41, 42|, can be applied. In what follows we use the abbreviations 

as = as{f^l = e-^^^l}) , = as(/u^o = e-^«//f_o) (4-24) 

and suppress all references to x or the Mellin variable A^. It is understood that in x-space all 
products of quantities depending on x have to be read as Mellin convolutions ([4.3[). 



The general solution of Eq. ( [4.2[ ) can be decomposed as 

= qLoM) + <lLM) (4-25) 

with the boundary conditions 

qLoM,o) = , qlMfi) = Q"(/^/,o) • (4.26) 
The homogeneous solution can be written as 

qIM) = ^(«s) f^)"""" U~\a,) q^ifil,) (4.27) 

with 

Ro = , U{a) = l + aUi+a^U2 + ... ■ (4.28) 







^We have taken the opportunity to correct the two misprints (one sign and one superscript in the coefticient) 
in Eq. (14) of ref. These misprints have also been noticed in ref. p^ . 
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The expansion coefficients Uk are complicated functions of the sphtting-function matrices i'C-'^) 
in Eq. ( ^ ) and the coefficients /3i<k of the /3-function ( p. 17 ) of QCD. A recursive exphcit repre- 
sentation (see ref. for a derivation) is given by 



k 



k 



(4.29) 



Here r± denote the eigenvalues of Rq in Eq. ( 4.2g| ) and e± the corresponding projectors, viz 



Rn 



r_e_ + r+ej 



(4.30) 



with 



1 

Wo 



qq 



gg 



3(0) 

gg 



(0) d{0) 



4PA.^P 



1 



-Rn 



rzpl 



r± - rzp 

The matrices R^ in Eq. (4.29) read Ri = Ri and 

k-l 

Rk>i = -Rfc + ^ Rk-iUi , Rk>i 



1=1 



/3o 



(3 



(4.31) 



(4.32) 



where ' represents the coefficients of a^~^^ in the second part of Eq. 



The coefficients Ui and U2 in Eq. ( [4. 281) are required at NNLO. The quantities Uk>2, on the 
other hand, receive contributions from beyond-NNLO splitting functions, hence these terms can be 
left out. It is understood that then also U~^{aQ) is expanded, and that all terms beyond second 
order in the coupling constant are removed in Eq. (4.27). In this manner the spurious poles in 
Eq. ( 4.29| ) at A^- values where r_ — r-|_ it vanishes are completely cancelled in the solution ( 4.27| ). 

The inhomogeneous solution gj^ijom 'with the boundary condition ( [4.26| ) is given by 



—1 7 / 

'^em tinhorn V*^! 



with 



= -Uias)a-^'' r daa^°~'^U~^{a)R^{a) 
k=o ^0 k=l V^o J 



(4.33) 
(4.34) 



where P^ ^ stands for the coefficients of Oem^s the first part of Eq. ( |4.5| ). Using Eqs. ( 4.2g| ) and 
( 4.34 ), the expansion of Eq. ( [4.33 ) leads to 



^em 9inhom(.A''/J 



l-(^ 



[i-r^)-^r:^ 



(l + as[/i)|l- ^"|i?o^(fi7-C/ii?o^) 



(4.35) 



l-(^ 



l--Ro" 



:i + R,)-^ [r] - U^Rl - - [/?] i^o") + . . . , 



where we have again only written out those terms which can be consistently [^3[ included at NNLO. 
The non-singlet solutions can be obtained from Eqs. ( [4. 27] ) and ( 4.35| ) be replacing all quantities 
by their scalar non-singlet analogues. In particular Eq. ( |4.29| ) is replaced by Uk^na = —\Rk,na- 
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5 NNLO quantities in x-space 



We now proceed to the explicit x-space results for the photonic coefficient functions and the photon- 
parton splitting functions required for the NNLO evolution. As far as our results are exact, we 
shall write them in terms of the harmonic polylogarithms -f^mi,...,m„ (a^), n^-j = 0, ±1, introduced 
in ref. [44| to which the reader is referred for a detailed discussion. The lowest-weight {w = 1) 
functions H^ix) are given by 

Hq{x) = Inx , H±i{x) = =Fln(l=F2;) . (5.1) 

The higher-weight (w > 2) functions are recursively defined as 

— - In^ X , if nil, ''^w = 0, . . . , 

(5-2) 

dz f.ra^{z) Hm2,...,m^{z) , else 

with 

Mx) = - , f±i{x) = . (5.3) 

X 1^ X 



In analogy to the notation for harmonic sums [45|, a useful short-hand notation is 

^^0,...,0,±i,0,...,0,±i,...(^) = ^±(m+i),±{n+i),...(a;) • (5.4) 



For w < 3 the harmonic polylogarithms can be expressed in terms of standard polylogarithms [46|; 
a complete list can be found in appendix A of ref. [^]. A Fortran program for the functions up 
to weight four has recently been presented 



For the convenience of the reader we recall, in this notation, the MS splitting functions and 
coefficient functions obtained previously |^, ^, ^, Henceforth suppressing the argument 'x' of 
the harmonic polylogarithms, the first and second-order splitting functions in Eq. ( [4.5| ) read 

(5;ipW(x) = 4pqg(x) , p^^ix) = X^ + {l-xf , (5.5) 

and 

5-^P^}^\x) = CF|56 - 116x + 80x2-4(l-4x)Fo-8(l-2x)Fo,o-16i/i (5.6) 
+ 16pc^g{x) (^-C2 + Ho + Hi + Ho,o + i^i,o + Hi^i + ^2) } 

CipW(x) = C7^|^- 64 + 32x + yx2-(24 + 40x)//o-16(l + x)Fo,o| , (5.7) 
where p = ns, q. The corresponding first-order coefficient functions in Eqs. ( [4.12| ) and ( p. 30 ) are 

C'cg'(^) = 5n.'cg"(^) = -4pqg(x)(4 + Fo + Fi) 

C'ci'i'(x) = 5ns'ci',^r(^) = 16^(1-^) • (5-8) 

All these results can be directly read off (by replacing Ca and 1/2 Nf — > 1) from their hadronic 
counterparts 4!g(x) fH, ||, Pig^\x) and P^^ (x) ||, in the latter case additionally keeping in 
mind that the off-diagonal quantity Pg^\x) does not contain a 6{1 — x) contribution ^. 
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The same holds for the NNLO photonic coefficient functions not written down before. These 
functions can thus be inferred from the results of refs. |0|, yielding 



Cf 



16 1294 478 72 



+ 32 1 X + -x^ + -x^ C2 

V 6 2 h r 



9 



+ 64 1 + -x^ Cs - -T - + 



16 /1 59 339 



15 \x 



+ — I ^ + 90 + 40x + 36x^ )i/_i,o-6( 1 



X + 81 I i^o - 28 I 1 



20 12 



7 



X + y X^ ) Hi 



16/1 



44 ,2 32 o , 
— x + 24x2 + yx3)//o,o 



40 36 o 

52 1 X + — x^ 

13 13 



(Hifl + - 32 ("l - ^x + ^x^] 



+ 32x |^C2 H^i + Hq + Hi + 2//_i,_i,o — — 4-^^0,0,0 — Hiflfl 

^H2,o - \h2,i - H3) - 32pqg(-x)(C2 H.i + 2F_i,_i,o - if-1,0 



— ~H2 n — 
2 ' 


^^2,1 - 


+ 32pqg(x 


) (C2 i^O 


3f. 

- 4^2,0 - 


- -f^2,l — 






ll5x 


256 



1 

4^ 



(5.9) 



^(x) 



64 /1 13 39 



Cf { X H x^ H X H x'^ C2 H \ x - 9x^ Hq 

ll5x 15 5 5 3 V 5 15V2; 4 2 / " 

16 (1 + 3x - 4x2) i^i +^ (^-^ - 5x + 6x3^ H.i^q ~ ^ ("^ + ^'^^) ~ 32x//2| 



(5.10) 



Analogous to Eqs. (3.2)-(3.7) of ref. Q and Eqs. (3.3)-(3.6) of ref. [|14[ for the hadronic second- 
order coefficient functions we also provide compact approximate representations of these results 
in terms of logarithms, 

Lo = Inx , Li = In(l-x) . (5.11) 
With deviations up to a few permille, Eqs. ( ^ ) and ( |5.10| ) can be parametrized by (p = ns, s) 

- {(26.67 - 317.5 (1 - x))L? - 72.00 - (1287x-^ - 908.6)Li + 1598 LiLq 
- 13.86L|] - 27.74Lg - 67.33 Lq - 1576 - 2727x + 3823x2 - 0.59(5(l-x)} , 

- {(-187.2 + 193.7x)(l - x)Ll + 737.0 (1 - x)Li + 452.2 LiLq 



84.87xLg - 64.25 (1 - x)Lo - 66.44 (1 - x)} 



(5.12) 



In the MS scheme the small 5(1— x) contribution in the first relation (of course absent in the exact 
expression, but useful for obtaining high-accuracy moments and convolutions) is not relevant in 
x-space at the present accuracy in the electromagnetic coupling. 
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The 0{a^^a'^) NNLO photon-parton splitting functions Pp^\x), p = ns,ps, g — which are not 
related to their hadronic analogues by any simple substitutions — are not completely fixed by our 
results in Sect. 3. Therefore we resort, for the time being, to approximations analogous to those 



derived in refs. [13, 14, 15 1 for the hadronic parton-parton splitting functions. In the MS scheme 
adopted for our calculations, the leading-A'f term of the NNLO photon-quark splitting function, 
for example, can be cast in the form 

^nsPl^ili^) = = E + /smooth (X) + X:^ni0 • (5-13) 

m=l n=l 

Here Li and Lq are the end-point logarithms defined in Eq. ( 5.1l| ), and /smooth is finite for < a; < 1. 
This regular function constitutes the mathematically complicated part of Eq. ( ^.13|) . See Eqs. (21)- 

and ( |5.10| ) above) to extract 



(25) of ref. |44] for a systematic procedure (not applied in Eqs. (|57 
the end-point terms of the harmonic polylogarithms entering the exact expression. 

At NLO the leading large- rr soft-gluon contribution, ln^(l — x), to Pnl\ cancels in the physical 



kernel ( 4.22 ). Anticipating the same cancellation at NNLO, the coefficient A4 in Eq. ( ^.13 ) can be 
inferred from the known splitting functions and coefficient functions. Besides this term we keep 
two of the remaining three large-x logarithms and two of the small-x logarithms in Eq. ( 5.13| ) , and 
choose a two-parameter ansatz (a low-order polynomial in x) for /smooth- These parameters and 
the coefficients of the selected logarithms are then determined from the six coefficients of a^N^ on 



the right-hand-sides of Eqs. (3.1). Varying all these choices we obtain about 50 approximations. 

The two representatives spanning the error band for most of the x-range, marked by 'A' and 'B' 

(2) 

below, are finally selected as our estimates for P^^i^ i and its residual uncertainty. 

Analogous procedures are applied for the {N^) pure-singlet photon-quark splitting function and 
the term of the NNLO photon-gluon splitting function. The non-singlet and gluon pieces 
are smaller in absolute size and uncertainty than the terms, hence for them it suffices to select 
just one central representative. The resulting approximations are displayed graphically in Fig. 1 
and Fig. 2. For the non-singlet case, the selected parametrizations (shown as solid lines) read 



128/27 Lj + 3.8636 + 97.512 Lj 



511.199 x + 84.808 Ln 



128/27 



22.878 Ll 



1319.749 X 

1 



+ 



5.4658 Lf - 295.331 Li - 1803.989 x 
+ 740.532 - 259.036 Lg + 27.110 L^} + 



^ns -fns7,2(^) 



(5.14) 



with 



5-lp(2) 



ns7,2 ^ 



73.672 



= Nf {-0.2949 Lf + 34.854 Li + 157.995 

- 33.059 Lo + 2.887 Lg} . (5.15) 
The corresponding approximations chosen for the pure-singlet splitting functions are given by 
S'^P^ll^j^ix) = iVf {(1-x) (-110.167x2 + 876.629 x + 23.605 X" 

+ 668.725 xLo + 387.125 xLg + 121.403 Lq} 



^p(^) (x) 



-1 



iVf {(1 - x) (-34.853 X + 305.244 + 50.950 x 

+ 101.246 xLq + 220.083 Lq - 5.0738 L^} . 



(5.16) 
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Finally the NNLO photon-gluon splitting function and its present uncertainty are parametrized by 

5~'^Pgl]A{x) = (1 -x) (769.616 Li + 1329.961 x2 - 391.569 2; 

+ 317.267 Lg + 60.519 Lg + 15.018 x^^Lo + 6^^P^^\{x) 



^s^Pgj,Bi^) = (1 - ^) (-105.632 Lf - 415.549 - 429.907 



357.604 Lo - 146.286 Lj] + 64.666 x'^Lq + 6~^P^l\{x) (5.17) 



with 



^s^P^^M^) = iVf 1(1- (43.748 Li + 61.028x2 - 70.910 - 33.914 



X 



- 105.172 Lo + 1.3972 L^} . (5.18) 

In all cases the averages 1/2 [A + B] represent our central results. 

Denoting the coefficients of of, on the right-hand-side of Eq. ( 4.20 ) by A^jg , the additional 



NLO (/ = 1) contributions to the DIS-y splitting functions are given by 

A^^P (x) = 4Cf {7 - lOx + (1 - 16x + 32x2) Hq + (g - 12x + 16x2) {H2 - C2) 

+ (5 - 36x + 32x2) ^ (2 _ 4^ + 8^2^ jj^^^ ^ 4pqg(x)(i7i,o + 2Fi,i^ 

(^^"i Aj]}^^ (x) = 8/3 Cf {2x-^ - 20 2x -M6 x2 - (3 -M5 X - 4x2) 

- (3 + 4x-^ - 3x - 4x2) ^ ^ ^^(^^ _ ^^^^ _ ^^)| (5 -^9) 

with p = ns, s. The corresponding NNLO terms can be parametrized as 

(5p ^ Agi^^P (x) = 9.482 Lj + (33.37 + 2585 (1 - x) ) Lf 122.6 + (5598 + 7949 (1 - x) ) Li 

- 9825 LqLi - 2.963 7.407 L|] - (176.0 - 2616 x) - 828.6 Lq - 1851 30120 x 

- 7595x2 - 0.65 (5(1 -x) 

iVf { - (1.044 + 32.57 (1 - x) ) Lf -h 26.75 - (0.266 - 615.1 (1 - x) ) Li -h 4.557 LqLi 

- 0.529 + (12.23 - 38.59 x) Lj + 41.91 Lq + 75.74 + 733.7 x - 1003x2 + 0.05 (5(1 -x)} 

+ iVf/ps {(2.083 - 68.96 Li)(l - x)^ - 86.40 (1 - xfloLi + 1.778 + 3.278 (5.20) 
+ (41.86 + 105.6x)(l - x)Ll + 1.241 (1 - x)2Lo - (15.80 x"^ + 94.39 - 8.281 x)(l - x)^} 

and 

C^AgJ<;^(x) = -{25.77 Ll + 766.7 Li)(l - x) - 1337 Lq^i - 4.741 + 4.741 

+ (83.11 - 464.9 x)Lg 443.1 Lo - (60.36x-^ - 1772 - 650.3x)(l - x) 
-FiVf { - (0.737 L? - 216Li)(l - x) 310.6 LqLi - (40.99 - 36.66x)Lg 

- 113.7Lo + (15.80x"^ - 331.3 - 48.88x)(l - x)| , (5.21) 



where fps has been defined in Eq. (p.2|). Eqs. ( 5.20| ) and (5.21) deviate by a few permille or less 



from the (somewhat lengthy) exact expressions deferred to Appendix B. 
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Figure 1: Approximations of the N^^ parts Pp^\(x) of the NNLO photon-quark (p = q) and photon- 
gluon (p = g) sphtting functions, as obtained from the six moments calculated in Sect. 3. The 
selected representatives (5.14) and ( 5.17| ) are shown by the full curves. 
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Figure 2: As Fig. 1, but for the Nl contributions Pp^2(.^)^ with p = ns, ps (left) and p = g (right) 



The full curves represent the parametrizations in Eqs. ( |5.15| ), (|5.16 ) and (|5.18| ). 
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6 Numerical results in £c-space 



In this section we finally present the numerical impact of the NNLO corrections on the evolution 
of the photon's parton distributions and on the structure function F2 ■ Concerning the evolution 
kernels we confine ourselves to the respective inhomogeneous contributions P'^ and I2 to Eqs. (p 



and ( 4.21 ). The higher-order homogeneous (hadronic) non-singlet and singlet kernels have been 



discussed in detail in refs. |T^, and refs. [14|, respectively; see also ref. [^] for a recent brief 
summary using the updated NNLO singlet splitting functions [15|. Also our subsequent illustrations 
of the solution of the evolution equations are restricted to the photon-specific inhomogeneous piece 
( 4.35 ) and its contribution to F2 ■ For brevity the results are shown only with all scales identified, 



i.e., at fir = fJ-f for the parton evolution and at fj,r {= fXf) = Q for F2 ■ The experimentally elusive 
longitudinal structure function will not be addressed here. 

In Figs. 3 and 4 the singlet photon-quark splitting functions Pq^ are displayed in the MS and the 
DIS^ factorization schemes. The NLO and NNLO curves are shown for as = 0.2 which corresponds 
to a scale between about 20 GeV^ and 50 GeV^, depending on the precise value of as{M'^). After 
removing the charge factors Jns and 6s defined in Eq. ( 2.291) only the NNLO terms depend on the 



number of flavours; the curves presented refer to A'^f = 4. In the MS scheme (Fig. 3) the NNLO 
corrections are small except for small and very large values of x, amounting to less than 2% for 
0.1 < X < 0.95. Larger corrections for x ^ 1 are obvious since the NNLO splitting function 
ln^(l — x) is more singular in this limit than its NLO analogue, P^^^ ~ ln^(l — x). Likewise 
large NNLO effects are expected for small x due to the first non- vanishing pure-singlet term ~ 1/x. 

The higher-order corrections to Pq^ are somewhat larger in the DIS^ scheme (Fig. 4) due to 
the absorption of the large coefficient function (72,^ (which, in particular, reverses the sign of the 
leading large-x contribution). Here the NNLO corrections, under the conditions specified above, 
reach +6% at x ~ 0.6 and exceed —6% for x > 0.9. In the right parts of both figures the relative 
NNLO effects are also displayed (dotted curves) for the non-singlet splitting functions Pns7, thus 
the pure-singlet effect can be directly read off from these figures. Also this effect is larger in the 
DIS^ scheme where it exceeds 1% at x < 0.3, instead of only at x < 0.15 in the MS scheme. 

The corresponding results for the photon-gluon splitting functions xPg-y are shown in Fig. 5. 
The relative NNLO corrections (not shown separately for these quantities) are larger than in the 
photon-quark cases. However, the absolute size of Pg^ is much smaller than that of Pq^ except at 
small X. In the MS scheme Pg^ remains negative at x > 0.2 at NNLO. In the DIS^ scheme, on the 
other hand, Pg^ is positive at large x, but seems to turn negative at NNLO for x < 0.05. 

Fig. 6 depicts, again for as = 0.2 and A^f = 4, the inhomogeneous contribution I2 to the physical 



non-singlet evolution kernel ( 4.21| ). Here the NNLO corrections are particularly small, about 1% 



or less for 0.05 < x < 0.95. As in Figs. 3-5, the present uncertainties arising from the residual 
error band for the O^a^^a^) photon-parton splitting functions are estimated by the NNLOa and 
NNLOb curves which derive from the upper and lower approximations, respectively, in Eqs. ( 5.14| ), 



( 5.16 ) and ( |5.17] ) together with Eqs. ( 5.15 ) and ( 5.18| ). These uncertainties are virtually negligible 



at X > 0.25 and remain perfectly tolerable for x > 0.05, where they amount to less than ±1% for 
Pq^ and I2 with respect to the central results ^ (NNLOa + NNLOb) not shown in the figures. This 
accuracy rapidly deteriorates towards small values of x. 

In Fig. 7 and Fig. 8 we present the inhomogeneous solutions ce^Jni ^/inhom' / = ^) ff) the 
singlet-quark and gluon distributions. The NLO approximation derived already in ref. § is ob- 
tained from the NNLO expression ( [4.35 ) by removing the third line, the Ui term in the second line 



and the U2 contribution to the first line; for the LO result also the rest of the second line and the 
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Figure 3: The perturbative expansion ( [4 .Bp of the photon-quark splitting function Pq^(x,as) in 
the MS scheme for Og = 0.2 and A'^f = 4. The relative NNLO corrections in the right part are 
also shown for the non-singlet splitting function Pns7- Here and in the following three figures the 
subscripts 'A' and 'B' indicate the respective approximations in Eqs. ( ^.14 ), ( 5.16| ) and ( |5.18| ). 
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Figure 4: As Fig. 3, but for the photon-quark splitting functions ( [4.20 ) in the DIS-.y scheme 
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Figure 5: The NLO and NNLO approximations ( fl.5| ) for the photon-gluon sphtting function 
Pg-y(x, Os) in the MS scheme (left) and the DIS-.,, scheme (right) for = 0.2 and A'f = 4. Note that 
x-P is displayed here, unlike in Fig. 3 and Fig. 4. 
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Figure 6: The perturbative expansion (f4.22| ) of the inhomogeneous physical kernel /2^(x, Cg) for the 



evolution (4.21) of the structure function ^2^^ = 1/x -^2^8 'default' value of the scale [ir 
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Ui term in the first line have to be ignored. Like the NLO illustrations in refs. 0, the results 
are shown at iJj = 50 GeV^ for |Ujq = 1 GeV^ and A^f = 3. For the strong coupling constants we 
use the realistic values as(?7ijg) = 47rag = 0.45 at LO and 0.42 at NNLO and NNLO, and solve 
Eqs. ( |2.16| ) directly at D = 4 with the appropriate number of terms in Eq. ( |2.17| ), i.e., including the 
coefficients Pk<m |38| for the N'^LO evolution. Here and in Fig. 9 the (barely visible) differences 
of the NNLOa and NNLOb curves include the present uncertainties (S^ of the hadronic NNLO 



splitting functions. The pattern of the NNLO corrections for the solutions roughly follows that of 
the corresponding splitting functions in Figs. 3-6. At this order fij = 50 GeV'^ is about the lowest 
scale where the inhomogeneous MS gluon density is positive over the full x-range. 

Finally the structure functions -^2^inhom(^' '3^ = 50 GeV^) resulting from these parton densities 
(and their non-singlet analogues) are shown in Fig. 9 for the MS and the DIS^ schemes. In both cases 
we have taken care to avoid spurious higher-order contributions which would arise from a simple 
convolution of Eq. ( [4.35| ) with the corresponding expansion of the hadronic coefficient functions, 
see Fig. 3 of ref. [^0| for a NLO illustration. Note that the boundary conditions for i^2''inhom 
are different in the two schemes: in the MS scheme this quantity is given by the corresponding 
photonic coefficient function ( [4.12 ) at = /^jq' "^here it vanishes in the DIS^ scheme according 
to Eq. ( [4.18| ). Thus, besides a 'physical' input describing the F2' at = ^^q, a large additional 
'technical' contribution to q'~'{fi'jQ) in Eq. ( 4.26| ) is required in the MS case. As mentioned above 
Eq. ( 4.21 ), the complete structure functions evolve approximately like non-singlet quantities at 
large x. In fact, under the conditions of Fig. 9, the complete DIS^ result for -F'2^inhom ^^^^ 
non-singlet approximation differ by more than about 1% only at x < 0.2, 0.25 and 0.3 at LO, NLO 
and NNLO, respectively. 
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Figure 7: The inhomogeneous LO, NLO and NNLO contributions to the photon's singlet-quark 
and gluon distributions in the MS scheme at /ij = 50 GeV^, as obtained from Eq. ( 4.35 ) for A'f = 3 
and /i^o = 1 GeV^ with asinj^) = 0.45 at LO and 0.42 at NLO and NNLO. 
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Figure 8: As Fig. 7, but for the inhomogeneous solutions in the DIS-y factorization scheme ( 4.2(]| ) 
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Figure 9: The perturbative expansion of the structure function F:^^^^^^ at = 50 GeV^, as 
derived from the results shown in the previous two figures and their non-singlet counterparts by 
means of Eq. ( [4.10| ) in the MS scheme (left part) and Eq. (4.18) in the DIS^ scheme (right part). 
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7 Conclusion 



We have calculated the next-to-next-to-leading order QCD corrections to electron-photon DIS 
and to the evolution of the photon's quark and gluon distributions. Our exact results for the 
corresponding photon-parton splitting functions are presently confined to the first six even-integer 
moments. Thus the practical applicability of the NNLO evolution is, for the time being, restricted to 
not too small values of the Bjorken variable x. This restriction seems to be more serious here than in 
lepton-hadron DIS, as the photon-parton splitting functions enter the evolution equations directly, 
not via smoothening convolutions with non-perturbative initial distributions. Consequently the 
effect of these functions is accurately known over a considerably smaller range of x than that of their 
hadronic counterparts; a residual uncertainty of about ±1% or less is found for the total photon- 
quark splitting functions at NNLO only at x > 0.05. However, only the perturbative component of 
the photon's parton densities is affected by this uncertainty, and while this component dominates 
at large x, it represents only a small correction to the hadronic (homogeneous) contribution at 
small X. We thus expect that our results are sufficient for extending NLO analyses like those of 
refs. [51-54] to NNLO for the full x-range covered by the measurements at LEP. 

Our present calculations are limited to effectively massless quarks, hence they do not apply to 
the charm and bottom contributions to the structure functions at small and intermediate scales. 
These contributions have been computed at NLO in refs. |55, corrected figures for F^^^^^^ 



have been presented in ref. [^]. The NLO effects are found to be fairly small for this quantity, 
indicating that the NNLO corrections may be rather negligible [^] . While a full massive NNLO cal- 
culation does not seem feasible at present, these corrections could be estimated using the threshold 
resummation as done for the lepton-nucleon case in ref. [^]. 

By confining ourselves to the lowest order in the electromagnetic coupling we have assumed, as 
usual also in QCD analyses of lepton-hadron DIS, that the QED radiative corrections are treated 
elsewhere. The corresponding formalism has been set up in ref. |^] for measurements of the 
photon structure via the process e^e~ e~^e~ + X. The non-factorizable corrections due to 
photon exchange between to two electron lines have been shown to be negligible in ref. |6^, using 
a pseudoscalar as a simple model for the final state X. The dominant corrections arising from 
photon emissions of the 'tagged' electron line have been investigated for realistic measurements of 
the photon structure in refs. Neither of these studies has addressed the QED corrections to the 
subprocess 7*7— We expect these corrections to be rather large, especially at the high scales 
which, hopefully, will become accessible to precise measurements in 67 collisions at the future linear 
collider. Such contributions can be included in our present formalism by extending the analysis to 
higher order of electromagnetism. This extension is left to a future publication. 

Fortran subroutines of our NNLO coefficient functions and of the approximations of the 



NNLO splitting functions can be found at http : / / www . nikhef . nl/ '-^avogt 
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Appendix A 



Here we present the analytic expressions for the anomalous dimensions and coefficient functions up 

to order a^^^ at the cvcn-intcgcr values N = 2, . . . , 12. The notation is as in Sect. 3; in addition 
Ca and Cf arc the standard QCD colour factors, Ca = Nc = 3 and Cp = {N^ ~ l)/i2Nc) = 4/3, 
and Q stands for Riemann's (^-function. The photon-quark anomalous dimensions are given by 
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The photon-gluon anomalous dimensions read 
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The photonic coefficient functions for the structure function F2 at fir = t-i-f = Q are given by 
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piv=io _ 8704 / 72533010722807 216 \ 

(-^paemj C2,:, - 7425 +««^^(^ 3486607740000 + 55 

9. ^ / 209277550337208434702179 4286102504597 

+ aUusCpNi 

s^ps ^ t 1^ 287046836661816000000 1981027125 



35264 
1981027125 11 



V zo/u^uooouuioiuuuuuuu lyoiuz/izu 11 

2 /_ 9568565675783568703229 366130385968 8856 200 ^ 

+ F 46259614172772000000 ^ 3679050375 ~ 3025^^ ^ "33"^V 

0^0/ 26284376777719892724358177 3565449892804, 17712, 2096, 

' V 117545679613013652000000 77260057875 3025 33 ^ 

2^ /462262990271149761829 302180374 , \ 

^ V 12721393897512300000 66891825 ^ J 
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._i piv=i2 _ 8110049 f 6258789011950819 44 

(^aem) C2,V - "7567560 ^ ''^ V 299266727760000 + 13^^ 

2 / 80246164545568080582270787397 _ 3109586855303 499488 \ 

+ fls/ps F f 62342050024259309304000000 847767375 ^ 91 

2 / 109517569276335628551465332167 59504875891511 2563 580^ 

+ F A I 517918569432308108064000000 615479114250 '"^ ~ 1014 ''^ 273''^ 



194382513719914205037949320835327 542538591728921 2563 4440 \ 

777654732002610624258096000000 ^ 16617936084750 ^ "SOf^^ ~ "gT^V 
2/1937462197313198708954658827 294410659 \ 

^ ^ \ 51791856943230810806400000 ~ 79053975 ' ^ 

and the corresponding results for Fl read 

a ^-1 P.^=2 232 32 \ 

(Jpaem) c^;^ = 3 ( ~Y35 ~ 

2^ ^ /^2686 3808^ 896^ 
2^ ^ / 11788 6512^ ^^^^ 

-^^-^-^(^-^^3) (A.19) 
5528^ 
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2 s r. ( 48227851 4384^ 

. / 51480649 174284 

+ "^^^^^ r^5437r " ^72^^^ + 128C5 

2^2 / 912595079 451816^ ^^^^ 
' V 39375 105^ V ^ ^ 



.-1 p,iv=6 2 / 137761\ 



2, „ / 97510126879 61778 1760 
+ aJpsCf ATf 4084101000 + ^iT^' - 

2„ „ / 43056771637 17446^ 480^ 
+ a^CpCA \- 437582250 " + — 

2 2 101548438631 87148 _ 960 ^ 
+ ""^^ I 1021025250 + ^20^^3 " — Cs 

2^ /826872491 64 \ , ^ 



31 



_1 p,Ar=8 o . _ / 51097\ 



[dpaem) = —+asCF 



2 / 166525045327877 617672 1664^ 

+ aJpsC^iVf 1620304560000 + 1^^' " 

2 / 103132581580249 393671 128 

^ V 1485279180000 ~ 779625 ~3~''^ 

/ 2384408424295187 15446822 ^ 256, 

' ^ V 35646700320000 779625 3 

2^ /36203618923 1216 \ , 

+ a^CpNf Cs A.22 

'V 5304568500 1485^ / ^ ^ 



._i p,7v=io _ 4 / 509195549\ 

(<)paeW c^,^ - - + asCp r 352182600 j 



2 / 551115536481640147 694547002 10240 

+ Os Jps F f 2761393330080000 ^ 1091475 '"^ 11~''' 
2 / 321273125948555419 290422414 320 

^ ^ V 6230251067040000 ^ 156080925 '"'^ ^ Tl'"^ 

2^2^9053411269935853949 1721766952 , 640, 

+ aiCi \ 

' ^ V 188465094777960000 156080925 ^ 11 

2^ 195749323625753 272 \ , 

+ alCpNi C3 A.23) 

'V 38073756520800 429 



197186301194 \ 



^•"^ "91 V~ 179847793125/ 

2 / 4911117022004232715153043 463564050788 126208 " 

+ Os/ps F f 15569942563501326000000 496621125 ''^ 91 ''^ 

2 / 47756006182490429534639 5631384334 1920 

'^^^ P \ 1197687889500102000000 2152024875 91 

2^2 / 281711596115081884853551 4695483724 3840 

' ^ V 7784971281750663000000 717341625 91 

2^ 780777857066540089 3392 \ 

+ atCpNi Cs • (A.24) 

' V 194429852191575000 6825^/ ^ ^ 



Appendix B 



The exact expression for the ©(Ogj^Og) contribution ( ^.20 ) to the transformation ( [4.20|) of the 
photon-quark spHtting function to the DIS^ scheme is given by 



2 1 o 72 



/431 4 16 102 

\ \ T. .1 



4 1 308 48 
48 + — ^ + —X + 76x2 ^ 3 j II _^ - 4(9 + lOx + \^x^)R-2,- 
15 3 5 ' 



1,0 
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/ , 2 1 130 „ n 24 o\ , /1883 4 1 1683 373 o\ „ 

- 2 (7 - 14x + 6^2) iJi,iC2 - 4 (3 - 6x + 7a;2) /fg^i _ 2 (s - 2x + Gx^) i72C2 

/ 9\ / 9\ 1679 4 1 4549 276 o 

-2(5-10x-4x>,C3-4(2-2x + 9x^)i/4-33^-^--^x + — 

/809 4 1 15521 213 o 32 o\ /5 2\ 

- (12O + 15 X + + + j ^"-'^ - (2 + + j 

- 2 (3 - 5x + 2x^) H2fi - 2 (3 - 6x + lOx^) iJi,3 - 2 (3 - 6x - 2x2) iJi,oC2 

y + 17x - 60x2^ H2,i - + + 16x2 + ^x^^ ^3 - Q + liar - 15x2^ Hi,o 

3 - 8x + 4x2) ili,2 - 2 (1 - 2x + 6x2) Hi^^^q^q - (l - 2x + 12x2) ^2,0,0 

1 287 „^ o 48 o\ „ /I „ „ 2\ /6 4 66 2\ .9 

2 + — + 36x2 ^ \ ^^^^^^ _ / _ + 2x - 8x2 j ^^^^^ \ ^2 

^ + 73x + 16x2 + ^x^^ + (x ~ + ^°^^) ^^'1 + + 36x2)//o,oC2 
+ 4 (3 + 2x + 8x2) ^ + 2 (1 - 2x - 2x2) iJi,o,o,o + (l - 2x + 16x2) iJ2,i,i 
+ 2 (3 + 6x + 4x2) (2/f_i,o,o,o + 2/^-1,3 + 8//_i,_iC2 + 8i7_i,_i,_i,o - 4i?_i,_2,o 

- 4F_i,_i,2 - 4i7_i,oC2 - 7//_iC3 - 6i7_i,_i,o,o) - 2 (l9 + 14x + 32x2) i?_2C2 

+ 8 (2 + X + 8x2) ^^^^ + 4 (^5 + 2X + 8x2) ^ + 4 (^7 + + 12x2) i7_2,0,0 

308 _ 2 -iA > /431 8 1 4328 102 2 32 o\ ^ 
+ ^13 + _. + 62x2 + 24x3jc3+(^ + ^- + ^x--.^ + Yx3jc2 

+ {^- 74x + 60x2^ + (^40 + + ^x + 44x2 ^ ^^3^ 

/ 4 1 92 48 o\ / 172 , 48 o\ 

+ (32 + —-2 + ya; + 12x2 + —x^\ H_i^2 + f 16 + — x + 76x2 ^ _^3 \ 

+ (Si/l,!,!,! - 2ffi,i,o - 6//l,l,l,o - 10//i,i,2 - 12ifl,2,0 - 14//l,2,l - 6i/2,l,o) Pqg{x) 

/613 8 1 4 1 4274 288 2 32 ,\ „ 1 
+ (1^ + 45^ + 15 X + + + 

^ ^ r / 11 1 394 ^ „ 44 o\ 233 111 12283 7 2 

/167 464 464 o\ /571 542 248 n 44 o\ 

701 466 433 ^\ /56 158 / 

— - —X + — x2 j + 22x + — x2 j //_2,o - 2 (2 + 5x + 3x2) 

95 191 433 n\ „ /40 122 88 o\ „ , (. ^ „ 2\ 

y - -^a; + -Y"^ ) ^2 - ( Y + + y ^ ) ^-1,0,0 - 4 (^1 + 2x + 3x2) 

/ . i\ , „ . „ ..X /1283 111 21301 3814 ^\ 
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/II 85 . / 1 3 22 o\ ,2 /II 121 29 n\ 

+ (3 - 6x + 8x2) ^^^^ ^ ^ ^ ^g^2^ ^ ^ ;L2x + 16x2) F_2,-l,0 

/II 19 o\ /34 107 79 n\ ^ /46 137 97 n\ , 

+ + y X + 6x2 j ^3 + _ + j ^^^^ + (- + -X + -x2 j H_,C2 

/95 1361 433 n 44 o\ , /68 214 158 o\ H 2.. 

+ y - + — x2 - -x3 C2 + y + — X + — x^ F_,_,o - yX^i^.i 



+ yH^i(s + 4f/'_i^_2,o — 8-H'_i^_iC2 — 8i?-i, -1,-1,0 + 6//-1, -1,0,0 + 4//-i,-i,2 + 4/7-1, 0C2 
-2F_i,o,o,o - 2i7-i,3) Pqg(-x) + 4x (/7o,o,o,o + ^4 - 2Fo,oC2) + f ^^iCs + 4Fi,oC2 



—2//i, 0,0,0 + 4-ffi,iC2 — 2/71,1,0,0 — 2/71,3 — 2/72,o,o — ^-f^i,i,o — 3" ^2,0 — ^^1,1,1 ) Pqg(x; 



yA^irC^ |(1 - 2x) {H0C2 - H3) + (3 - 6x + 2x2) i7i,o,o - Q 
+ (4 + 8x + 4x2) (77-1,0,0 - ^-iC2 - 277-i,-i,o) - (4 - 8x + 4x2) 77iC2 



- - X + 2x" ) 77o,o,o 



/41 44 , 80 2 24 3\ ^ / 55 70 2 24 3 



+ X H x^ x-^ 77o - 7 X H x^ x"* C2 

V12 3 3 5 y ' V 3 3 5 



83 2 1 1091 553 2\ rj f^^ 80 80 2 



20 15x 60 15 / V 3 3 3 



X + —^x^ 77o + -X + — x^ 77 



'1,0 



,„ 70 2\ 83 2 1 7 69 2 /55 76 70 2 » 

+ 7 - 13x H X'' 772 H \ \ x x^ + x H x^ 77i 1 

' 3 7 10 15x 60 5 V6 3 3 ' ' 



+ (11 - 6x + 18x2) Ca + 2x2/72,1 + (277i,i,o + 77i,i,i + 2772,o) Pqg(x) + 8 (l + x2) 77-2,o 

2 1 16 24 o\ /137 105 125 

i H K -\ X H x'^ 77-1 + X H a 

15 x2 3 5 / ' V 12 2 3 



+ 16/psiVFCF |2(1 - 2x) (^774 - 77oC3 - 77o,oC2 + ^^o,o,o,o - ^Cl) + (l - 4x - y x2^ 77o,o,o 

175 20 1 4327 4132 . / 40 1 850 ^\ 

+ X x2 - 3 33x + x2 77i 

27 27x 27 27 V 27x 27 / 

+ (^1 + 71x + ^x2^ 77o + (^7 + 33x - ^x^^ 77o,o - (^4 + llx + ^x^^ (C2 - ^2) 

l + 8x-^x2Vc3 + 77oC2-^3)| , (B.l) 



where the function Pqg{x) has been defined in Eq. (5.5). The corresponding result for the NNLO 
transformation of the photon-gluon sphtting function, parametrized in Eq. ( 5.21| ), reads 

C^Ag)gS(x) = 

00^2 r/ 67 4 1 209 4 3\ ^ / 4 8 2\ rr ^36 26 \ 2 

/ 85 13 1 25 2 2\ /, 22 \ ^ ^ /I 29 , 2\ 

+ (-^ + y- + ^x + -x2 j 77i,i - (5 + -xj 77oC2 + + y^ " ^o.o.o 

/ 5 4 1 13 o\ ^ / 5 2 1 5 2 2\ 

+ X + 2x Ca + ^ h -X x^ 77i 

V2 3x 2 y^ V2 3x 2 3 / 

/ 9 326 2 2 4 3\ ^ /I 11 1 4 8 2\ ^ ^ 

+ X x2 + -x"* C2+ -H X x'^ 77i - (7 - x) HnG 

V 4 9 3 5 V3 3x 3 3 7 ' ^ ^ ^ 
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+ 4(1 - X) ( H^2C2 + 2H_2,-l,0 - H^2fi,0 + IhiC2 ) + (1 + x) ( 3F_iC2 - 3^^-1,0,0 



+ 6-f/"_i,_i,o + — 3i^o,oC2 — -f^2,0,0 — 2-^2,1,0 " 7-H'2,l,l — 2/^3,1 + --^0,0,0,0 + 3-f/'4^ 

12x2^ Hi+(^6 + 5x + ^x^^ i?2,o 

x3^ i/0,0 + (5 + ex)H3 + (^^ + + y X^^ //2,1 
+ (//-1C2 - ^-1,0,0 + 2/7_i,_i,o)pgT,(-x) + (^^1(2 - i^l,l,0 - P97(^) 



2 2\ /41 13 1 15 

+ l- + ^^^ + r)^^+(l2 + 12x + Y^ 

'39 2113 8 2 4 
+ 'T + ^^-3^ -5 



/3389 4 1 311 64 
I 1 X X 

V 360 45 X 40 5 / 



+ ( ^^^^ + + ^x - — x^ ) i?o + + — - - —X - — x^ - 8F- 



2179 149 1 77 39 

3,0 



180 180 X 15 5 



, ,'20 2 1 2 16 o\ /19 8 1 25 n , 

+ 32CaCf -hr+o-+oa;- — x^ H.q - i ^ - ^- + + Ax^ ] i/2,1 



3 3x3 3/ ' V3 3x 3 

11 14 4 2\ /31 2 1 128 2\ 
\ X + -x^ i^s \ X - 17x2 Hi 

3 3 3yv9 3x9 y 

8 4 1 5 8 2\ rr /37 67 35 2 



3 - 3 X + r + r) - U + 1^" + ¥^ j 

-, r 4 2\ /I 8 1 5 2 2\ /7 2 1 95 56 2 , rr 

1 - 5x x^] H^2 + h -X + -x^ 1 + X x^ H2 

3 / ' V3 3x 3 3 / ' V9 3x 9 9 ' 

~ (9 ~ ~ ''^ (-^-1-1,0 - -H'-iAO - ^^-1,2 + ^i?-iC2^ Pg7(-2^) 

+ (^/?i,o,o + 2Fi,i,o + 2i?i,2 + 3/7i,i,i - ^i/iC2) Pg7(a^) + (y + y ^ + ^^') ^oC2 

+ (1 — x) (^3f/'„2C2 - 2-?/_3,o + 2//_2,-i,o — 2i7_2,o,o — 2i7_2,2 " -f^oCs - 2//o,0,0,0 

+ (1 + X) (2/72,0,0 - 3F2C2 + 4/^2,1,0 + 6/^2,1,1 + 4/^2,2 + 4^3,0 + 6^3,1) 

/21 27 \ 2 /19 103 2\. /7 51 1 7 2\ 

+ \ X C? + \ X + 2x^ Ks + h -X x'^ 

VlO 10 / ^^ V6 6 y^'^ V2 3x 2 3 / ' 

/97 67 128 2\ rr 317 1 1 295 571 2 . ^ ^ 

+ ^ X x"' i^o H X H x'^ + 2x (Fo - Hq 0C2 + 2i74 

V18 9 9 / 18 6x 6 18 ^ ' 

32,, ^ r/11 26 1 1 10 ^\ 127 10 1 35 152 2 

H /V^F^F { h -X x'' Hi \ X H x^ 

3 lV3 9x 3 9 / 9 9x 9 9 

10 2\ rr /13 7 8 2A rr ^ / 16 22 4 2\.. 

3 + 21x + — x^ j Hq+^— + -x + -x^j (C2 - H2) - (^y + y ^ + 3^ j ^0,0 
+ (i7i,o + 2/7i,i)pg^(x) + 2(1 + x) {^Hs - C3 - ^oC2 + i^2,o + 2/^2,1) } ■ (B.2) 
Here we have used the abbreviation 

4 ^ 4 2 

Pgiix) = — + 1 - X- -X . 

For the numerical evaluation of the harmonic polylogarithms up to weight four entering Eqs. (|B 
and ( [B.2D we have employed the program of ref. |47|. 
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